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Chronopharmacoldnetics of nicotine 

Sjiti^Hale: For hlgh-clearilnce drugs such as nicotiiie, hemodynamic changes throughout the day may be 
expected to influence the rate of metabolism. 

Material afid tnethods: To assess the eflects of meats and diurnal rhythms on nicotine clearance^ an 
intravenous infusion of nicotine bitartrate was administered for 48 hours to II subjects. Two models 
to determine nicotine clearance variation throughout the day arc described. Both models were used to 
estimate the mean effect of meal and diurnal rhythms on nicotine clearance and individual parameters 
that were regressed against baseline covariates. Olearaucc was modeled as a function of time {CL(t)] 
and split it in three components: a (constant) baseline value (Bi)} its circadian (diurnal) variation, and 
the effect of meal: CLi(t) = [Oj + circadiati(t)] [1 + meal(t)]. A two-compartmental (time-variant) 
model incorporating CL(t) was then fitted to the data providing estimates of CIi(t) conditional on 
literature values of the time-invariant parameters (volume of distribution and intcrcompartmental 
clearances). 

Raallst The estimated circadian(t) showed a maximuia at approximately 11 am and a flat minimum from. 
6 PM to 3 am; the estimated meal(t) showed a sharp increase up to 1 hour (after the meal), at which point 
clearance is increased 42%, and a slower decrease thereafter, returning to baseline (zero) after 2.8 hours. 
Individual estimates of baseUne clearance are found to have a linear relationship with body wei^t. No 
other covariate, sex in particular, effect could be found- (Clin Pharmacol Ther 1996;60:385-95.) 


Jean-Michel Grics, PhamiD, Neal Beaotritz, MD, and Davide Verotta, PhD 

San-Francisco, Calif. 


Smokers tend to regulate the intake of nicotine, 
adjusting their smoking behavior to maintain a cer¬ 
tain concentration of nicotine in the body.^ Con¬ 
versely, levels of nicotine in the body may influence 
when a person smokes a cigarette, how much nico¬ 
tine he or she inhales from the cigarette, or both. 
The level of nicotine in the body is determined by 
the dose taken in and the rate of elimination. Nic¬ 
otine is eiimiimted primarily by hepatic metabolism 
by way of C-osidation to cotininc,^ Nicotine is me¬ 
tabolized rapidly, so that the hepatic extraction in- 
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dicates that the clearance of nicotine is likely to be 
dependent on liver blood flow. 

Changes in drug metabolism throughout the day 
have been described for various drugs,^’* For faigh- 
clearance drugs such as nicotine, hemodynamic 
changes throughout the day may be expected to 
influence the rate of metabolism. For example, eat¬ 
ing a meat increases liver blood flow and may in¬ 
crease drug clearance. We have shown previously 
that eating a meal lowers blood nicotine concentra¬ 
tions/ consistent with an increase in clearance. Car¬ 
diac output and liver blood flow vary throughout the 
day as a function of activity with a decline during 
sleep-^ It might be expected that nicotine clearance 
would decline overnight as a function of reduced 
liver blood flow. Consistent with this hypothesis is 
the observation that during constant intravenous 
infusion of nicotine for 24 hours, there is a progres¬ 
sive increase in plasma nicotine concentrations 
throughout the night® despite half-life values after 
brief infusions of 2 to 3 hours. This observation 
suggests the possibility that nicotine clearance is 
reduced overnight compared with that during the 
day. 

This study aims to characterize the chronophar- 
macokinetics of nicotine In people. A 48-hour con- 
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tinuous intravenous infusion was administered to 
examine blood levels of nicotine throughout file day 
under steady-state dosing conditions. Resultant 
pharmacokinetic data were modeled to determine 
the effects of meals and circadian influences on the 
clearance of nicotine. 


Data analysis. Generally nicotine concentration 
kinetics are well described by a biexponential func¬ 
tion,® which, can be written (in differential form) as 
follows: 


dA.(t) 

dt 




(la) 


METHODS 

Subjects. Eleven healthy subjects (six men and five 
women) who were 23 to 54 years of age and who 
smoked 20 or more cigarettes per day were studied. 
They ranged in weight from 57.5 to 94.0 kg (mean 
70-4 kg) and in height from 155 to 182 cm (mean 
171,5 cm). Five (45%) subjects were white, four 
(36%) were black, and two (19%) were Asian. Two 
additional subjects dropped out of the study because 
they could not abstain from smoking for 48 hours. 
Baseline covariates (age, sex, weight, height, and 
race) were recorded for the pharmacokinetics anal¬ 
ysis. All subjects gave written consent after being 
advised of the nature of the risks of the study. Tlie 
subjects were confirmed to be in good health by 
physical examination, medical history, and clinical 
laboratory tests. Subjects were admitted to the Gen¬ 
eral Clinical Research Center at the San Francisco 
General Hospital for 5 days, during which time 
smoking was not permitted. Compliance with non¬ 
smoking was assessed by periodic measurements of 
expired carbon monoxide concentrations. Meals 
were administered at fixed times each day for all 
subjects. Breakfast was given at 8 am, lunch at 12 
noon, and dinner at 5 pm. Caffeine-containing bev¬ 
erages were prohibited. 

Experimental protocol. After a washout period of 3 
days, an intravenous infusion of nicotine bitartrate 
(1.25 mg base/hr) starting at 8 am was initiated and 
continued for 48 hours. One catheter on one arm 
was used for infusion, and another catheter on the 
other arm was used for blood sampling. Blood sam¬ 
ples were taken for measurement of nicotine con¬ 
centrations before infusion was performed, every 2 
hours for the next 22 hours (until 6 am), then every 
hour for the next 20 hours (until 2 am), and then 
eveiy 2 hours for the remaining 6 hours (ending at 
6 am). In. addition, a sample was taken 30 minutes 
after each mean was begim (i.e., 8:30 am, 12:30 pm, 
and 5:30 pm) to better characterize the kinetics re¬ 
lated to meals. 

Concentrations of nicotine were measured by gas 
liquid chromatography with nitrogen-phosphorus 
detection as described by Jacob et al.® and modified 
fox the use of a capillary column. 


dA,(t) 

dt 



Vs*-V, 


-A,(t) 


(lb) 


in which and are the amotmt of nicotine in the 
central and peripheral compartment, respectively, 
Vc is the volume of the central compartment, Vgs is 
the volume at steady state, Q is the inter- 
compartmental clearance [Q = kjiiVjs — Vc), 
where is the elimination constant from compart¬ 
ment 2 to compartment 1], and CL the constant 
clearance of nicotine from the central compartment. 
In this article we modeled clearance as a function of 
time (CL(t)) and split it in three components: a 
baseline vdue, the circadian (diurnal) variation, and 
the effect of meal. The model used for clearance can 
be written as follows; 

CL(t) = [hi + circadian(t)] [1 + ineal(t)] (2) 


in which represents the baseUne component, cir- 
cadian(t) the diurnal variation of clearance, and 
meal(t) the fraction of increase in clearance caused 
by the meal Substituting CL(t) for CL in equation 
la obtains a time-variant (versus time-invariant) 
two-compartment model, which will turn out to be 
more appropriate to describe nicotine kinetics. 

Mostly to compare across two possible ap¬ 
proaches, we developed two different models: a 
semiparametric and a parametric one. For the para¬ 
metric model: 


circadian(t) = ejsin 



(3) 


in which to, the period of the circadian component of 
clearance, was fixed to 24 hours, and Dj (amplitude) 
and 03 (phase) were estimated from the data. The 
effect of meal is given by the following: 


0604 


meaKt)= 


ift>tk 

otherwise 


(4) 


in which t^^ indicates the time of the k-th meal, and 
04 , 05 , and 06 are model-specific parameters (in this 
case a Bateman function). 

For the semiparametric model we used nonpaxa- 
metric functions (splines instead of the cosine and 
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exponential functions), but maintained a parametric 
twt5-compartmeiital structure. Splines are flexible 
functions that are basically piecewise polynomials 
which join at certain locations called breakpoints (£) 
and match their derivatives up to the degree minus 
one of the polynomials. For example, a linear spline 
matches only the values of the polynomials at the 
breakpoints to give the somewhat familiar “broken 
line.” We used cubic splines in this article. 

For the circadian component we used a periodic 
cubic spline; this matches all derivatives at bound¬ 
aries of the period (in this case 0 and 24 hours) in 
the same way as the cosine function used previously. 
The number of free parameters in a periodic spline 
equals the number of its breakpoints minus 2, For 
meal(t) we used the following: 

f2>(t-tj ift>h 

mcal(t) = < k (5) 

[ 0 otherwise 

in which s(t) indicates a so-called natural cubic 
spline (this has second derivatives at the first and 
last breakpoint equal to zero). To make s(t) a mean¬ 
ingful model for the meal effect, we constrained it to 
be positive and equal to zero at the time of meal and 
at the last breakpoint. The position (in time) of the 
last breakpoint of the spline is one of the parameters 
to be estimated; it characterizes the duration of the 
meal effect because s(t) reaches zero at this time and 
thereafter. The number of free parameters in the 
meal effect spline equals the number of breakpoints 
minus one. 

For both splines the first breakpoint was posi¬ 
tioned at zero, and the other breakpoints were equi- 
spaced between the first and the last ones (at 24 
hours and at the estimated end of meal effect for the 
periodic and natural spline, respectively). To deter¬ 
mine how many breakpoints to use (increasing 
breakpoints increases the number of parameters in 
the splines and therefore their flexibility) we tried all 
combinations from two to six breakpoints and de¬ 
termined the best combination with the Akaike in¬ 
formation criterion.^® (This criterion requires a de¬ 
crease of the objective function of at least two points 
per additional parameter for selecting the larger 
over the smaller model.) 

In a first stage we fit AiIV,. as given by equation 
(1) (with CL(t) replacing the constant CL) to the 
pooled data from ^ individuals to estimate circadi- 
an(t), meal(t), and 0i, fixing Vc, Vgs, and Q to 
literature values (see Benowitz et al.® and following 


text). This was done for the two time-vaiying clear¬ 
ance models, parametric and semiparametric. 

In a second stage to characterize interindividual 
variability we fit the following model to eacli indi¬ 
vidual’s data; 

circadianj(t) = vijCircSaian (t + (6) 

where j indicates the j-th individual, and circadian 
(t) is the estimate obtained at step 1. 11118 model 
states that individuals differ in the amplitude of 
circadian effect and in, for example, the time of 
inaxiinum effect that is changed by the time shift 
For the meal effect we used; 

meaJjft) = (7) 

where meal(t) is the estimate obtained at step 1, and 
? indicates the breakpoints of the meal spline. Mul¬ 
tiplying the breakpoints by has the effect of 
scaling the function on the time axis; this allows a 
different time of maximum effect and of end of 
effect Finaliy, individuals’ estimates of baseline 
clearance (-vsj) were allowed. Models of this sort 
(often referred to as “template models,” where, e.g., 
me^(t) is the template that can be scaled and 
shifted in time to characterize individuals) are often 
used in the literature.**'^ The total number of pa¬ 
rameters to be estimated for each individual is five. 

From ciicadian(t}, meal(t), and the individual es¬ 
timates of 7 jj, i = 1,... we computed the more 
intuitively meaningful parameters: (1) T„a(j, the 
time of the maximum effect of meal, (2) the 

maximum effect of meal, (3) the time, after 
Tn,„j, at which X% of maximum effect is present, (4) 
Amjxj, the maximum of circ^iar^(t), and (5) Aminj, 
the minimuni of circ^ianj(t), tthich we will report, 
together with the individual estimates of baseline 
clearance, to express interindividual differences- 

To determine whether the measured covariates 
(sex, race, age, weight, and height) influence these 
estimated parameters, we used a Wilcoxon test for 
sex, Kiuskal-Wallis for race, and linear regression 
for continuous covariates effect (weight, age, and 
height). 'The computer software NONMEM*^ was 
used to estimate the parameters. We used a propor¬ 
tional error model for all fits. 

RESULTS 

Fig. 1 shows the fit of the data to a two- 
compartment model with = 63.8 L, Vss = 212.6 
1, Q = 183.1 L/hr, and a constant clearance of 71.6 
L/hr, as reported by Benowitz et al.® 
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Ftg. L Two-compartment body model fit to data versus dock time. The solid line is the 
prediction. The dots are the observed concentrations for individual subjects. The dashed lines 
connecting the dots arc for the male subjects, and the dotted lines arc for the females. The 
vertical long dashed lines are the meal starting times. 


Semiparametric model The fit of the semiparamet- 
ric model is shown in Fig. 2. There is a substantial 
improvement of the fit to the data compared with 
the two-compartment model with constant clear¬ 
ance (shown in Fig. 1) as assessed either by the 
objective function (which drops from 1682 to 1533) 
or by visual inspection. (Note that a drop of 14 
points would be sufficient to prefer the semipara- 
metric model.) The circadian spiine requires four 
coefficients (six breakpoints) and showed a maxi¬ 
mum at approximately 11 am and a flat minimum 
from 6 PM to 3 am (Fig. 3, A, solid line). Similar 
results were obtained with three and five coeffi¬ 
cients, although in both cases the flat region is sub¬ 
stituted by a single minimum (not shown). The 
spline describing the meal effect requires three co¬ 
efficients. The meal effect lasted for 2.8 hours with a 
sharp increase up to 1.1 hour and a slower decrease 
thereafter. Fig. 3, B shows the pooled estimate of 
meal effect (solid line) and the individuals’ estimates 


(dotted lines). The maximum effect of a meal was 
estimated to be a 42% increase in clearance. Fig, 3. 
C shows the pooled estimate of clearance versus 
time over the 2 days of the experiment (solid line) 
and the envelope including the prediction of all 
individuals (dashed line). (We do not show individ¬ 
uals’ estimates because the graph is too cluttered.) 
The effect of meal seems to be the most important 
component of clearance, either when we look at the 
plot of clearance versus time or when we consider 
the drops in objective function when we include one 
component of clearance at a time. 

Individual estimates of the parameters 

Tss,!. and are reported in 

Table I. The coefficient of variation of the baseline 
clearance is approximately 28%, similar to what has 
been previously reported.® The other parameters 
exhibit a higher variability: coefficient of variation of 
35% for and T,,,,, 45% for and A,„„, and 
66% for The mean estimates corresponding 
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Fig. 2. Semipararaetric fit of data. Sec leg&nd to Fig. 1. 


to the meal effect axe physiologically reasonable; the 
maximum increase in clearance is obtained 1 hour 
after the beginning of a meal, and the meal effect 
lasts 2.7 hours. Fig. 4 shows the predictions for each 
subject with their estimated parameters. Note the 
overall good fit to the data (which might be made 
better with a more flexible individuals’ model) and 
the high variability in the shape of the individual 
responses. 

Parametric model. The fit of the parametric time- 
vaiying clearance model shows the same (signifi¬ 
cant) improvement compared with the (time- 
invariant) two-corapartment fit. 

Individual estimates of the parameters Tn,jx,j> 
(which CQuaiS — An^gxj), 
and 75 j are reported in Table II. The main charac¬ 
teristics of the parametric estimate of the circadian 
component are similar to those of the semiparamet- 
ric estimate; the peak clearance is obtained at ap¬ 
proximately 10 AM, and the maximum decrease in 
the circadian clearance are of the same absolute 
value (7.4 L/hr for the parametric model versus 6.4 
IVhr for the semiparametric). The major differences 


between the two estimates of circadian component 
are that (1) the parametric estimate is symmetric 
around the baseline clearance, which is not the case 
for the semipararaetric one, and (2) the parametric 
estimate does (can) not show the flat portion of the 
circadian component seen in the selected semipara¬ 
metric estimate. 

The parametric estimate of the meal effect differs 
to a greater extent from the semiparametric esti¬ 
mate; the estimated clearance function is very sharp, 
with a maximum at 10 minutes and a slower decrease. 
The effect of meal at 214 hours after the meal is, 
however, almost identical; 4.1% and 3.5% for the 
semiparametric and parametric models, respectively. 

Computing the individual parameters with the 
parametric model leads to somewhat different re¬ 
sults (Table II). The most important differences are 
for the mean baseline clearance (68 L/hr for the 
parametric model versus 76 L/hr for the semipara¬ 
metric) and the effect of meal (69.6% maximum 
increase in clearance for the parametric model ver¬ 
sus 51.1% for the semiparametric). Sometimes the 
parametric model estimates T„ax to the time 
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Fig. 3. A, Estimated circadian component of clearance [cjrcadian(t)] versus clock time. The 
solid line is the pooled data estimate, and the dolled lines are the individuals’ curves. B, Estimated 
meal effect on clearance [meal(t)] versus time (in hours). Meal effect is expressed as a fractional 
inaeasc. The solid line is the pooled data estimate, the dotted lines the individuals’ curves. C, 
Estiinated resulting total clearance [CL(t)] versus clock time. The solid line is the pooled data 
estimate, the dashed-line is the envelope including all individuals’ predictions. The vertical long 
dashed lines arc the meal starting times. Estimates are based on the semiparametric model. 


of meal (particularly in subjects 2 and 9). The reason 
for this apparent lack of practical identifiability of 
(compare with the estimates obtained by the 
semiparametric model) is probably due to the sam¬ 


pling design (the first observation is half an hour 
after the meal). The estimates of T^j^ indicate a 
too-large peisistence of meal effect (the mean is 9.4b 
hours), which is one of the reasons to reject the 
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TabJ« I. Individual parametej: estimates (semiparametiic model) 


Subject No. 

A 

T 

•* maxj 

(hr) 


Tyyj 

(hr) 

(Uhr) 

A 

•^maxj 

(Uhr) 

1 

65.6* 

0.735 

0.421 

2.06 

-4.81 

4.53 

2 

49.2 

0.452 

0.391 

1.26 

-8.6 

8.1 

3 

81.7 

0.753 

0318 

2.11 

-736 

7.12 

5 

99.9 

0.483 

0.185 

135 

-5.02 

4.72 

6 

593 

1.44 

0359 

4.02 

-8.63 

8.33 

r 

593 

1.19 

1.42 

3.34 

-Z41 

2.27 

9 

102 

139 

0.284 

3.90 

-4.36 

4.x 

10 

54^ 

0.977* 

0.438* 

Z73’ 

-333 

332 

It 

103 

1.07 

0308 

2.98 

-10.9 

10.3 

12 

63.1 

0.852 

0.625 

238 

-7.47* 

7.04* 

13 

98.1 

1.27 

0367 

336 

-12.1 

1L4 

Mean 

76.0 

0.965 

0.511 

Z70 

-6.85 

6.45 

SD 

213 

0342 

0.336 

0357 

3.09 

2.91 


■Y 5 J, The baseline individual clearance; the time from beginning of meal to the maximal po$tcncal change in dcaiance, the onaximutii change 

in cleaisace cau»4the effects of the me the time, after Tnut^ at which X% of the individual meal effect is left; Anraj^ the. maxunum circadian 

decrease in clearance below the baseline clearanoe; Ag,^^, the maximum arcadian mcrease in clearance above the baseline clearance. 

*Median value for each parameter. 


parametric model in favor of the semiparametric 
one. 

In terms of goodness of fit the semiparametric 
and the parametric model give similar results. The 
objective function for the pooled fits is 1533.4 for the 
semiparametric model and 1534.5 for the parametric 
model. For the individuals’ fits the square root of the 
mean squared differences between individuals’ data 
and predictions is 0.106 for the semiparametric model 
and 0.107 for the parametric model. 

Coyariaies ^ecis. Individual estimates of baseline 
clearance are found to have a linear relationship 
with body weight (slope = 1.369; p = 0.006). No 
other covariates were found to be related with the 
estimated parameters. In particular, there is no sex- 
related difference either on the diurnal variation in 
clearance or in the effect of meals. 

DISCUSSION 

Our study is novel in two aspects. First, we pro¬ 
vide new information on the human biologic char¬ 
acteristics of nicotine. Second, we present an inno¬ 
vative modeling technique for studying meal and 
circadian influences on drug clearance. 

We have previously reported that eating a meal 
while receiving a continuous intravenous infusion of 
nicotine results in a lowering of piasma nicotine 
concentrations.® This observation is consistent with 
increased clearance caused by an increase in hepatic 
blood flow, which is known to occur after eating.^'^ 


In our earlier work we reported a change in nicotine 
clearance but made no attempt to quantitate the 
magnitude of change and duration of change in 
clearance. The observation of a progressive increase 
in plasma nicotine concentrations overnight during 
a 24-hour intravenous infusion of nicotine suggested 
a decline in nicotine clearance overnight, but the 
change in clearance was not quantitated.® 

This study was designed specifically to quantitate 
the temporal changes in nicotine clearance through¬ 
out the day with a particular emphasis on changes 
that occur with meals and at night. A 48-hour infu¬ 
sion was used to distinguish the accumulation of 
nicotine overnight caused by a long terminal half- 
life versus an actual change in clearance. Meals were 
administered at fixed times of day, and nicotine 
blood concentrations were sampled frequently after 
meals to be able to characterize the details of the 
meal effect on nicotine kinetics. 

The modeling techniques we use are to the best of 
our knowledge unique. We not only use nonpara- 
metric functions to represent diurnal variation, but 
we also incorporate the additional feature of mod¬ 
eling the effects of discrete events (meals), which 
appears to be previously unreported in the litera¬ 
ture. More standard modeling of circadian variation 
is generally limited to the use of cosine functions. 

The use of splines with different numbers of 
breakpoints as nonparametric functions and the 
Akaike criterion to choose the number of break- 
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Fig. 4. Individuals’ data and predictions of plasma concentration of nicotine versus clock time 
(individuals are ordered left to right, top to bottom). The dots are the observations, and the 
solid lines are the predictions obtained by semiparametric model. 
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points is also an important and useful impro\'ement. 
In this context it allov.'s us to test the complexity of 
the circadian and meal effect components, fewer 
breakpoints being associated with lower completdtyj 
without assuming a priori functional forms for these 
components. 

The results obtained by the semiparametiic and 
parametric approaches are similar with respect to 
the size of the meal effect and the esitent of circadian 
variation. The semiparametric estimate indicates 
that the circadian variation is possibly more compiest 
than a simple cosine function, obtaining a nonsymmet- 
ric function diowing a Hat region in clearance between 
6 PM and 3 am. The spline function used to represent 
the meal effect represents a slightly more realistic 
model than the parametric one, because it estiniates a 
discrete end of the meal effect instead of an a^mptotic 
attenuation. Individual estimates of the duration of 
meal effect are unrealistically long (see Table 11). 
The parametric estimate shows a somewhat unrealistio 
sharpness in the meal effect, with a maxiraum reached 
at 10 minutes after the start of the meal However, the 
sampling rate immediately after a meal is at most evciy 
half hour, which is insuffidenl to readt a definitive 
conclusion on the details of the shape of the meal 
effect. 

In general we prefer the simplicity of representa¬ 
tion offered by spline functions, their clean hierar¬ 
chic framework (in terms of increasing number of 
breakpoints/parameters), and what we consider a 
vital Safeguard in avoiding model misspecification 
(spline functions do not include any a priori knowl¬ 
edge and thus are influenced only by the data). Of 
course more complex (and realistic) parametric 
models can be used to represent clearance variation 
(for example, sums of cosine functions or sums of 
exponentials in equation 3), and such “parametric” 
models (we use quotes to stress their lack of physi¬ 
ological interpretabLity) can ofGer simitar advantages. 

We found that the average clearance of nicotine 
in this study over a 48-hour period was similar to 
that previously reported after a SO-minute infusion 
of nicotine.^ The effects of meals on nicotine clear¬ 
ance as illustrated in Fig. 3 were substantial, with an 
average 42% increase in clearance peaking at ap¬ 
proximately 1 liour after the beginning of the meal. 
The peak of clearance estimated by the parametric 
model occurs much sooner, and it is difficult to 
decide which estimate is the better one. Based on 
the semiparametric estimate the effect of meals 
lasted for approximately 3 hours. As a high- 
clcarance drug the clearance of nicotine is expected 
to be influenced by liver blood flow, and this is 


Table IL Individua] parameter estimates 
(parametric model) 


Subjea 

No. 

(L/i-) 

(hr} 

^thhXj 

(hr} 


1 

64.0'' 

0.553 

0.329 

4.13 

-4.65 

2 

40.6 

0.050 

0-671* 

11.7 

-9J9 

3 

65.8 

0.503 

0.710 

15.2 

-12.8 

5 

97.5 

0.251 

0.320 

1.96 

. -3.12 

6 

43.9 

1.31 

0.852 

17.4 

-10.0 

7 

45.3 

0.251* 

2.41 

10.9 

-10.1 

9 

102 

0.050 

0.244 

9.42* 

-5.88* 

10 

53.2 

0.603 

0.347 

5.45 

-3.61 

If 

92.4 

0.050 

0377 

12J 

-14.2 

12 

60.0 

0.402 

0.721 

7.66 

-5.11 

13 

87.2 

0.955 

0.715 

7.33 

-2,41 

Mean 

63.0 

0.452 

0.700 

9.46 

-7.41 

SD 

22.7 

0.398 

0.606 

4.72 

4.09 


* Median value for each parameter^ 


presumed to be the basis for the effect of meals on 
the clearance of nicotine. Eating a meal produces an 
average 40% increase in iivet blood flow,’'*^ similar in 
magnitude to the change in clearance of nicotine. 
The metabolism of other drugs has been shown to 
be increased by meals, although in no case has an 
estimate been made of the magnitude or duration of 
the change in drug clearance. 

There was a significant circadian variation in 
nicotine clearance, but the magnitude of the cir¬ 
cadian change in clearance (approximately 17% 
from peak to trough) was much smaller than the 
effect of meals on nicotine clearance. The times of 
day when peak and trough clearances were ob¬ 
served are of interest. The peak was seen at ap¬ 
proximately ll:00 AM and the trough from 6 pm to 
3 AM. As discussed previously the clearance of 
nicotine is likely to be influenced by liver blood 
flow. Circadian changes in liver blood flow mea¬ 
sured as the clearance of indocyanine green have 
been reported.^ Liver blood flow is highest at 8 am 
and lowest in the late evening and early morning 
hours. Liver blood flow is also increased bj' physical 
activity, which is greater in the daytime compared with 
nighttime hours. Thus chaises in liver blood flow can, 
explain the circadian pattern of increased daytinae and 
decreased nighttime clearance of nicotine and other 
high-clearance dnigs such as midazolam,^ bupiva- 
caine,'* and teibutaline.^® 

Also of note in our study was the considerable 
individual variability in the extent of meal-related 
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and circadian changes in plasma nicotine levels (and 
therefore clearance) as shown in Fig. 4. Some sub¬ 
jects had quite substantial fluctuations in nicotine 
in association with meals and substantial day/night 
variation, whereas others had relatively flat 
curves. The reasons for individual variability in 
meal and circadian effects in nicotine clearance 
are not clear but might be related to individual 
differences in sieep/wake cycle, level of physical 
activity, or both. 

Our findings have several potential biologic im¬ 
plications. Circadian changes in nicotine clearance 
could influence cigarette smoking behavior. Smok¬ 
ers typically smoke after meals, and a decline in 
plasma nicotine concentration caused by an increase 
in clearance could contribute to the urge to smoke. 
Some smokers smoke less in the evenings, and that 
behavior might be contributed to by the lower clear¬ 
ance in the evening compared with the morning. It 
was noted that there is considerable individual vari¬ 
ability in the extent of circadian variation in nicotine 
clearance. Conceivably this could contribute to in¬ 
dividual variability in temporal patterns of cigarette 
smoking. 

Circadian changes in nicotine kinetics have impli¬ 
cations in understanding blood nicotine concentra¬ 
tions during nicotine medication therapy, particu¬ 
larly with transdcrmal nicotine delivery ^sterns. 
Pharmacokinetic analyses of nicotine absorption 
riom nicotine patches have assumed a constant 
clearance of nicotine over a 24-hour period.® Be¬ 
cause nicotine clearance is actually higher in the 
daytime and lower in the evening compared with the 
mean, the estimated release rates for nicotine from 
the patches are actuaily higher during the day and 
lower at night than originally reported. 

Finally, that the clearance of nicotine is slower at 
night helps explain why significant plasma nicotine 
concentrations are seen in many smokers when they 
awake in the morning, even if they have not smoked 
during the night.Persistence of nicotine levels in 
the body overnight may have implications in under¬ 
standing tolerance to nicotine effects and potential 
nicotine-related toxicities. 

in summary, we present a novel analytic approach 
for characterizing changes in drug clearance 
throughout the day. We find significant meal-related 
and circadiaD changes in nicotine clearance, which 
could be important in understanding the human 
pharmacologic characteristics of tobacco use and 
nicotine-containing medications. 
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